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Abstract 



■ We examine the nonequilibrium dynamics of a self-interacting A^ 4 scalar field 

{Sj ■ theory. Using a real time formulation of finite temperature field theory we 



derive, up to two loops and (3(A 2 ), the effective equation of motion describing 
the approach to equilibrium. We present a detailed analysis of the approxi- 
mations used in order to obtain a Langevin-like equation of motion, in which 



the noise and dissipation terms associated with quantum fluctuations obey a 



fluctuation-dissipation relation. We show that, in general, the noise is colored 
(time-dependent) and multiplicative (couples nonlinear ly to the field), even 
though it is still Gaussian distributed. The noise becomes white in the infi- 
nite temperature limit. We also address the effect of couplings to other fields, 
which we assume play the role of the thermal bath, in the effective equation 
of motion for (j). In particular, we obtain the fluctuation and noise terms due 
to a quadratic coupling to another scalar field. 
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I. INTRODUCTION 



The possibility that the Universe went through a series of phase transitions as it expanded 
and cooled from times close to the Planck scale has been actively investigated for the past 
fifteen years or so It is hoped that by studying the non-trivial dynamics typical of the 
approach to equilibrium in complex systems, many of the current questions of cosmology, 
from the origin of the baryonic matter excess to the large-scale structure of the Universe, 
will be answered in the near future. As is well-known, the origin of density perturbations 
that seed structure formation has been linked to either the existence of topological defects, 
such as strings or textures formed during a GUT-scale transition 0, or to inflation in one 
of its incarnations. In particular, the old, new, extended, and natural models of inflation 
all invoke a symmetry breaking transition in which nonequilibrium conditions play a crucial 
role || . At the electroweak scale, the focus has been in generating the baryon number excess 
during a first-order phase transition [|J. Even though there are certain questions related to 
the reliability of the perturbative expansion for weak enough transitions || as well as to the 
mechanism by which weak first-order transitions complete 0, it is currently believed that 
nonequilibrium conditions are a crucial ingredient for baryogenesis. 

Despite its relevance, not much has been done to understand nonequilibrium aspects 
of phase transitions in cosmology. (This situation is rapidly changing. We will soon refer 
to past and recent work on the subject.) Most of what has been done so far is related to 
the finite temperature effective potential (computed in general to one-loop order) which, 
by its very definition, is only adequate to describe equilibrium situations; the calculation 
is usually done in euclidean time so that we can obtain the equilibrium partition function 
from a transition amplitude. The great advantage of using the effective potential is that 
it gives us information about static properties of the system such as its possible stable 
and metastable equilibrium states, and critical temperatures for phase transitions. The 
disadvantage is that we lose all information about real-time processes, which are crucial to 
understand the mechanism by which the system approaches equilibrium. In fact, the one- 
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loop effective potential does carry, in a somewhat indirect way, information about unstable 
states in the system. These are states which are in the "spinodal" region, where the effective 
potential is concave. If we start with the system in thermal equilibrium above the critical 
temperature and then quench it to below the critical temperature so that its order parameter 
takes a value within the spinodal, the approach to equilibrium will be initially dominated 
by the growth of small amplitude long wavelength fluctuations, in the mechanism known 
as spinodal decomposition. Thus, the effective potential tells us that some states will be 
unstable, and that their final equilibrium state is at its global minimum, but it does not tell 
us how the system gets there. The reader is referred to the recent work of Boyanovsky, Lee, 
and Singh for details [[?]. 

These limitations of the effective potential were pointed out by Mazenko, Unruh, and 
Wald, in work where they argued that for strong enough couplings, the slow-roll approx- 
imation necessary for successful inflation may not be adequate. Instead, the approach to 
equilibrium would proceed by the formation and growth of domains, typical of spinodal 
decomposition ||. It was subsequently shown within the context of the new inflationary 
model, by both analytical and phenomenological numerical methods |Tt| that due to the 



small couplings needed for the generation of density fluctuations, the slow-roll picture of 
inflation was correct. 

This discussion of the validity of the slow-roll approximation in inflation raises some 
very interesting questions related to the way we picture the approach to equilibrium in field 
theories, which are quite independent of inflation. For example, the distinction between the 
"system", which is out of equilibrium, and the "thermal bath", which drives the system 
into equilibrium, is somewhat blurred in the context of nonlinear field theories. In fact, 
for self-interacting field theories, the short wavelength modes can serve as the thermal bath 
driving the longer wavelength modes, which have slower dynamics, into equilibrium. In this 
sense, the field can be its own thermal bath. Of course, other fields coupled to the order 
parameter scalar field (henceforth the "system") may serve as the thermal (or, at T = 0, 
quantum) bath. 



In one of the original works on this subject which was motivated by cosmology, Hosoya 
and Sakagami obtained an approximate dissipation term in the equation of motion satisfied 
by the thermal average of the scalar field, by invoking a small deviation from equilibrium in 
the Boltzmann equation for the number density operator. This calculation was then supple- 
mented by a computation of transport coefficients using Zubarev's method for nonequilib- 



rium statistical operators flll| . Using an approach which is closer to the one we will adopt 
here, Morikawa obtained the effective Langevin-like equation (that is, with both fluctuation 
and dissipation terms but not quite as simple as the Langevin equation) for a scalar field 
interacting with a fermionic bath using real-time field-theoretical techniques at zero and, 
very briefly, finite temperature l^j . More recently, Hu, Paz, and Zhang analyzed the case 



of a quantum bath given by a scalar field quadratically coupled to the system while 



Lee and Boyanovsky considered the case of a thermal bath given by a scalar field linearly 



coupled to the system [14]. Some works dealing with nonequilibrium evolution within a 



cosmological framework can be found in Ref. [HJ. Here we will only be concerned with 
dynamics in Minkowski spacetime. 



Recently, and in particular in Refs. |TB[ and |TJJ , the properties of the noise as being in 
general colored and multiplicative (unless the coupling between system and bath is linear) 
have been emphasized. This can have very important consequences to our understanding 
of phase transitions, as suggested by Habib, even though results at this point are prelim- 
inary ||16|| . The reason is that potentially, a multiplicative noise may sharply decrease the 
relaxation time-scales in the system and thus accelerate the approach to equilibrium. Nu- 
merical simulations of the approach to equilibrium have so far employed a phenomenological 
Langevin equation, with white and additive noise to mimic the effects of the thermal bath. 
In (l+l)-dimensions both the thermal nucleation of kink-antikink pairs ]I7[] and the decay 



of metastable states |TL8| were studied, while in (2+l)-dimensions the decay of metastable 



states was recently investigated |l9fl . The time-scales measured in these simulations agree 
with the theoretical prediction for the decay rate, T ~ exp[— B(T)/T], as long as B(T) 
is the classical (i.e., obtained with the classical potential) nucleation barrier given by the 



energy of the appropriate field configuration that saturates the path integral, the mass of 
the kink-antikink pair or the energy of the bounce configuration in the examples mentioned 
above. 

The question then is if the phenomenological Langevin equation used in the above simu- 
lations is indeed reproducing the essential physics of the approach to equilibrium, or if we are 
dangerously oversimplifying things. The above discussion suggests that the effective equa- 
tion which describes the approach to equilibrium of the slower moving modes can be quite 
different from the phenomenological Langevin equation with its white and additive noise. 
Two tasks are at hand then. First we must obtain the effective equation for a self-interacting 
scalar field which acts as its own bath and compare it with the equation obtained by having 
another field act as the bath. This should elucidate the nature of the thermal bath in these 
two situations, and also give us an answer as to whether the phenomenological Langevin 
equation is at all valid in some limit. The second task follows naturally the first. Once we 
have an effective equation we trust (in some limit), we should use it to simulate numerically 
the nonequilibrium dynamics, measure the relaxation time-scales, and compare the results 
with the results obtained with the simplified phenomenological Langevin equation. 

In this paper we will concentrate on the first task. Namely, we will obtain, within per- 
turbation theory, the effective equation of motion describing the approach to equilibrium of 
a self-coupled scalar field. We will integrate out the short wavelength modes whose influ- 
ence will be felt as a thermal bath through the nonlinear couplings to the longer wavelength 
modes, which we take as the system. The separation between bath and system is imple- 
mented by perturbation theory, since the effective action is obtained by integrating over 
small fluctuations about the state we are expanding about. We will include corrections up 
to two-loops, as nonvanishing viscosity (and transport coefficients, in general) terms in finite 
temperature field theory only show up by considering higher order (loop) corrections to the 
field propagators and are dependent on the imaginary part (decay width) of the self-energy 
corrections ]TT],^U|-|2I| . Fluctuation terms are obtained by associating the imaginary terms 



in the effective action as coming from the interaction of ip with fluctuating (noise) fields, as 



done in [0 and [14}|. We will also obtain the effective equation of motion in the presence 



of another scalar field quadratically coupled to the system, thus reproducing (even though 



we focus more on dynamical aspects) the analysis of Ref. (13| for finite temperatures. Our 
results could, for example, be used in the numerical investigation of symmetry restoration 
at finite temperature |2B] . 



The paper is organized as follows: In Sec. II we derive the effective action for a nonuni- 
form time-dependent background field configuration tp(x, t), up to two loops and order A 2 . In 
Sec. Ill we obtain the effective equation of motion for tp(x, t) and discuss the approximations 
involved in order that it obeys a Langevin-like equation. In Sec. IV we examine the effect of 
other fields interacting with the scalar field, by studying the case of a quadratically coupled 
scalar field and by evaluating its contributions to noise and dissipative terms. Conclusions 
are presented in Sec. V. Two appendices are included in order to obtain some technical 
results used in the paper. 

II. THE TWO-LOOP FINITE TEMPERATURE EFFECTIVE ACTION 

Consider the scalar field model with Lagrangian density 

m = ~ m 2 - - ^ (2.i) 

and with generating functional Z[J], in terms of an external source J, given by 

Z[J]= [D<f>e W {iS[<f>,J}} , (2.2) 

J c 

where the classical action is given by 

3[<f>, J]= f d A x {£[(/>] + J{x)<j>{x)} . (2.3) 



c 



In ( |2.3| ) the time integration is along a contour suitable for real-time evaluations, which we 
choose as being Schwinger's closed-time path |23,31],25], where the time path c goes from 



— oo to +oo and then back to — oo. The functional integration in ( |2.2| ) is over fields along 
this time contour. As with the euclidean time formulation, the scalar field is still periodic 



in time, but now with 4>(t,x) = 0(t — i/3,x). Temperature appears due to the boundary 
condition, but now time is explicitly present in the integration contour. 

As usual, the effective action T[(f] is defined in terms of the connected generating func- 
tional W[J] as 



W[J] - / d A x.J(x)f(x) 



(2.4) 



with <p(x) defined by f(x) 



«M arid 

5J{x) ' anQ 



W[J] = -iln / D<f>exp{iS[<f), J]} 



(2.5) 



The perturbative loop expansion for T[tp] is obtained by writing the scalar field as — > 
0o + f], where 0o is a field configuration which extremizes the classical action S[<p, J] and r\ 
is a small perturbation about this configuration. By using (|2.4|) and ( |2.5|) , we can relate 0o 
to if (0o = <P ~ v) an d write the effective action, for J — > 0, at one-loop order, as 



where 



r[tp] = s[<p] + -Trhi[n + v»(<p)] 



Trln [□ + V"((f)] = -i\n^Dr}ex-p{~ri[0 + V"((p))rj} 



(2.6) 



(2.7) 



Negleting contributions to (|2.6|) which are independent of f, we can expand the logarithm 



' _ ]Vm+l 



in as 



-Trln[D + n^)] = ^Tr c £ 

Z A m=l 



m v \ 2 



■ +oo (_-l\m+l / \\ rn n 

- £ — U Tr / • • • ^ 4 ^^ lA (xi - x 2 ) \<f 2 (x 2 )] G;» h (x 2 - x 3 ) . . . 



. . . \f 2 {x m )] G n r lm [x m - x x ) \f 2 (x 1 ) 



(2.8) 



The matrix representation in (|2.8|) is a consequence of the time contour, since now we must 
identify field variables with arguments on the positive or negative directional branches of the 
time path, that we denote by f + and <f-, respectively. As a consequence of this doubling 
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of field variables, we also have that G^\x — x'), the real-time free field propagators on the 
contour, are given by (l,n = +, — ) [^I] , ^6| 

i(T + 4>{x)4>{x')) 
i(T^(f)(x)(f)(x')) 

i{<P(x)<P(x')) , (2.9) 



G+ + (x-x') 
G^~{x - x') 

g;-(x-x') 

G^(x-x') 



where T + and T_ indicate chronological and anti-chronological ordering, respectively. G^ + is 
the usual physical (causal) propagator. The other three propagators come as a consequence 



of the time contour and are considered as auxiliary (unphysical) propagators p6 |. The 
explicit expressions for G^\x — x') in terms of its momentum space Fourier transforms are 
given by |2~T| , |2~5 



GAx - x) 



d 3 k 
(27)3* 



D ik.(x—x') 



G+ + (k,t-t>) G+~(k,t-t 
yG^+(k,t-t') Gj-(k,t-t>) ) 



(2.10) 



where 



G+ + (k, t-t')= G>(k, t - t')6(t - f) + G<(k, t - t')6(t' - t) 
G-f (k, t-t') = G> (k, t - t')6(t' -t) + G< (k, t-t')9(t- t') 
G+-(k,t-t') = G>(k,t-t') 
G^ + (k,t-t') = G<(k,t-t') 

and, for free propagators at finite temperature, 

1 



(2.11) 



G>(k,t-t')-- 



G<{k,t-t')=G>{k,t'-t) 



1 + 2n(w)) cos[cj(t - If)] - i sin[u;(t - £')] 



(2.12) 



where n(uj) = (e^ u — l) is the Bose distribution and lu = uj(k) is the free particle energy, 



u>(k) = yk 2 + m 2 . 
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Let us now add to Q2.6Q contributions up to two-loops and order A 2 . Graphically we 
have, 



™=*M + >o + >o< + >oo + -e- +o(as 



(2.13) 



where, in the graphic representation, (p is in the external legs and the internal propagators 
are given by Gr"? . In terms of the field variables ip + and <£>_, the terms in Eq. fl2.13|) are 
given by (note that now time runs only forward) 



%>] = J d 4 x {C[ip + ]-C[ip-\} 



(2.14) 



>0 



4 J J (2tt) 3 



Gt + (q,0) G+-(q,0) 



\ I (p 2 + (x) 



-<p 2 _(x] 



) 



<p+(x) - <fi(x) 



d 3 q 



[1 + 2n(u)] 



(2.15) 



(2tt) 3 2u(q) 

where G 1 !' (q, 0) is given by ( p. 11] ) (for t — t' = 0). Eq. (|2.15|) gives just the finite temperature 
mass contribution to the effective action (renormalized by a proper mass counterterm Sm 2 
which we are not including here). The second graph in (|2.13| ) is given by 



XX 



d 3 k 



i — Tr / d 4 xd 4 x' . 
16 J J (2tt) 3 



,ik.(x—x') 



d 3 q 



G+ + (q,t-t') G+-(q,t-t')) 
G7 + (q,t-f) G7-(q,t-t') 



x 



<fl(x' 

v 
A 2 



^ ( G+ + (q-k,t-t>) G+-(q-k,t-t>) \ ( 



i — / d 4 xd A x' 
16 



d 3 k 



G^ + (q-k,t-t') G7-(q-k,t-t') 



J 







-<pL[x) J 



jk.(x-x') I d 3 q 



?l{x)G+ + {q, t - t')G+ + (q - k, t - t') V l(x')- 



(2tt) 3 J (2vr) 3 

- V 2 + (x)Gi-(q,t- t')G+-(q- k,t- t'V-OO - V 2 -{x)Gl + {q,t- t')G7 + (q- k,t - 



+ v 2 _(x)Gr(q,t- t')G~(q- k,t- t')^ 2 _(x') 



(2.16) 



Equivalently, we get for the third graph in ( 2.13Q the expression 



XX) 



^(l + 2nM) [ dt , [J^M 



(2tt) 3 2u(k) 



(27T) 



G+ + (q,t-t'] 



G- 



t-f) + <p 2 _(x) GT-^t-f) 



(2.17) 



whereas the fourth graph in (2.13) is given by 



.-^ 2 f ,4 i4 i f d 3 k ik (x-x 1 ) f d 3 qi d?q 2 d 3 q% -> _ _ 
i — / d xd x / - — — e ( > / , \ , \ , \ S(k - q\ — q 2 - q 3 ) x 
12 7 J (2vr) 3 J (2tt) 3 (2tt) 3 (2tt) 3 V y y W 



x t - t')G+ + (q 2 , t - t')G+ + (q 3 , t - f>+0O" 

- t - t')G+~(q 2 , t - t')G+~(q 3 , t - t')^(x')- 

- <p_{x)Gj + (qi,t- t')G^ + (q 2 ,t- t')G^ + (q 3 ,t- t')<p + (x')+ 
+ ip4x)G^^,t-1/)G^ 2> t-1rjG^-{^t-1/)ip.(x , )} . 



(2.18) 



Before continuing, it is advantageous to rewrite the field variables tp + and </?_ in (|2.13|) 
in terms of new field variables ip c and (p&, defined by 



<P-= Vc - 



(2.19) 



The physical meaning of these variables is suggested in Ref. ||27|| , with cp& being basically 
associated with a response field while ip c is the physical field, which "feels" the fluctuations 
of the system. The change of variables ( 2.19Q will allow us to identify, in the effective 
action, the terms responsible for the fluctuations in the system (the imaginary terms). The 
association of ip c as the physical field imposes that we take (p& = (<p + = <pj) at the end 
of the calculation |21jl2||. In terms of the new variables ip c and if a, using (|2.11| ) and ( p,12|) , 
we get the following expression for the effective action (|2.13|) , using the physical propagator 
G^ + (q,t — f) in the Feynman diagrams ( |2.14| )-( pT8D , 
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., A r d 3 k (1 + 2n(u)) 



2 J (2tt) 3 2u(k) 



2 J J (2tt) 3 
A 



Gt + ^t-t')] 2 6{t-t') r d * k (l + 2n(-)) 



(2tt) 3 2cu(fc) 



V9 c (x)- 



- (4^ A (a;)^(x) + <p 3 A (x)(p e (xj) j + 
(i 3 A; 

+ 4<^ A (a;)^ c (a;)^(a; / ) 



4! 

+ f d*xd 4 x' 



ik.(x—x') 



d 3 q 



Im 



\J (2vr) 3 



A 2 

—fA(x)ip c (x')lm 

A 2 



V^A(a;)v5 c (a;)v?A( a;/ )+ 
G++(g,t - f - M - f)] *(* " 0~ 

0(t -0*(£ -ft 



d 3 g 



+ 



A 2 



Sik-qt-qz-qs)} . (2.20) 



LJ (2vr) 3 ' 

The last two terms in r[<^A<^c], Eq. ( |2.20|) , give the imaginary contributions to the effective 
action at the order of perturbation theory considered. It is straightforward to associate the 
imaginary terms in ( j2.20| ) as coming from functional integrations over Gaussian fluctuation 
fields 6 and £ 2 [0,0 

P>£iP[£i] / £>£ 2 P[£ 2 ]exp{z J d 4 x [<pa(x)<p c (xMx) + ^0*060*0 



exp | i / d xd x 

v2 



' ^2 

z— v?A(a;)v5 c (x)Re 



A 3 
+ «T77</?A(V)Re GJ + ,^a(0 



ipA(x')ip c (x') + 

(2.21) 



where P[£i] and P[£ 2 ], the probability distributions for £1 and £ 2 , respectively, are given by 



P[< L ] AT 1 <>xp { ™ / rf 4 xrf 4 x'6(x) ( y Re 



P[&] = iV- 1 exp { ~ J d A xd 4 x'£ 2 {x) (yRe 





2 




x,x' j 






G </> 


3 ) 

x,x' I 



6 GO 



(2.22) 



(2.23) 



where iVj 1 and A^ 2 1 are normalization factors, and in (|2.21 )-(2.23) we introduced the com- 
pact notation, 



11 



G 



d 3 k 



exp 



ik.ix — x' 



(2^r 



G++(q,t-t')G+ + (q-k,t-t') (2.24) 



.3 <P 



and 



Gl, 



_fk_ 

(2tt) 3 



exp 



- g x - <f 2 - ?3) 



^2.25) 



Therefore, using ( |2.21| ), Eq. ( |2.20| ) can be rewritten as 

T[<PA, ¥>c] = ~ In J £>£iP[£i] / £>£ 2 P[6] exp {^ cfr [y? A , y? c , £ 1; 6]} , (2.26) 

where 



Seffk^^b^] = Rer[^ A) ^ c ] + / d 4 a; [v?a(^)v ? c(^)^i (ac) + <^ A (^)6(^)] > 



^2.27) 



and Rer[<y9A, <p c ] is the real part of Eq. ( |2.20| ). In ( |2.27| ), the fields £i and £2, with probability 
distributions given by ( |2.22j) and ( |2.23| ), respectively, act as fluctuation sources for the scalar 
field configuration (p. £1 couples with both the response field ipA and with the physical field 
tp c , leading to a coupled (multiplicative) noise term (y>c£i) in the equation of motion for ip c , 
while £2 gives origin to an additive noise term. In the next section we examine the relevance 
of each of these noise terms in the equation of motion for the physical field (p c and evaluate 
the dissipation coefficients associated with them. 



III. THE EFFECTIVE EQUATION OF MOTION 



The equation of motion for tp c is defined by 

5ip A 



Iva=0 



. 



(3.1) 



Using JCTp and (^20|) , we obtain 



, A r d 3 k 1 + 2n(u) L , 

□ + m 2 + - / - — =M 1 - A 

2 J (2vr)3 2u(k) V 

A 2 r c l 
+— (f c (x) / d 3 x' / dt'Lp 2 c (x',t')lm 

Z J J —00 



4 *7A> 
-00 7 (2vr) 3 

A 2 



, + — d 3 x' / dt'(p c (x', f)lm 



(3.2) 
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where 



Gt 



and 



Gj 



3 

x,x' 



are given by (12.241) and (p. 251), respectively. In order to 
t' ' 

obtain a Langevin-like equation, a series of approximations must be performed in the above 
equation of motion. These approximations will certainly limit the scope of applicability of 
the final equation to be obtained (very much as in linear response theory), but on the other 
hand will elucidate important aspects of the nonequilibrium physics. Strictly speaking, a 
Langevin-like equation can only be used to describe the nonequilibrium dynamics of slowly 
varying modes in near-equilibrium situations. To see this, we now focus on the last two 
terms on the left hand side of Eq. ( p.2|) . 



A. Dissipation Coefficients 



Let us first consider the term in the equation of motion dependent on Gt + . In- 

L ™ J x,x' 

specting (|2.24 ), it is clear that the spatial nonlocality can be handled by considering only 
contributions with zero external momentum, as in the computation of linear response func- 



tions | 22"|| . This is what is usually done in the computation of the one-loop effective potential 
as an expansion of vertex functions with zero external momentum, which is physically equiv- 
alent to considering only nearly spatially homogeneous fields. We thus obtain, 



dt'(pl(x',t')lm 
A 2 



Gj 



2 
A 2 



dt' 



dt' . 

oo J (2?r) 3 



fl(x,t') - V l(x,t) 
d 3 q 



d 3 q 



Im 



G+ + (q,t-t') 



Im 



Gt + (q,t-t') 



+ 
(3-3) 



where we have summed and subtracted in ( |3.3| ) the last term in the rhs. In order to handle 
the temporal nonlocality let us further assume that ip c varies sufficiently slowly in time, so 
that we can expand the first term in the rhs of ( |3.3| ) to first order around t. This is a valid 
assumption for systems near equilibrium, when ip c is not expected to change considerably 
with time. (This has been called the quasiadiabatic approximation in Refs. |T2J and |]28|j.) 
We then obtain, 
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A 2 



<p e (x) / dV / dt'cpl(x',t')lia \GV 



A (p c (x,t)(p c (x,t) 



A 2 



dt' 



dt'(t' - t) 
d 3 q 



d 3 q 
(2^f 



Im 



(2tt> 



-Im 



Gt+(«f,t-t') 



+ 

(3.4) 



The emergence of a time direction within this approximation is surely related to ne- 
glecting the faster moving modes in the description of the dynamics. This is an interesting 
question which deserves further study, but that we will not address in the present work. The 
last term in the left hand side of Eq. ( |3.2|) can also be worked out as in ( ft. 3D and (|3.4j) and 
we obtain 



A 2 r '•' 
y / dV 



A 2 



dt'<f c (x',t')lm [GJ 
b c (x,t) f_Jt\t'-t)lm f[ / -^G^fat-t' 



A 2 



+—<p c (x, t) 



dtfhn 



3 

n 



(2tt)* 

G+ + (q v t-t>) 



(2vr) 3 * KHJ 



S(qi + <? 2 + 9*3) + 
5(9i + 92 + 93)- (3-5) 



The first term in the rhs of (|3.4j) and ( ft.5p are the corresponding dissipative terms 
associated with the fluctuation fields £1 and £2, respectively. The last term in the rhs of 
( ft.4j) is the one-loop finite temperature correction to the vertex [second graph in (|2.13| )1, 
while the last term in the rhs of ( ft.5|) is the contribution to the finite temperature two-loop 
correction to the mass coming from the "setting sun" diagram [the last graph in (|2.13|) 1. 
The time integrations in (|3.4j) and ( ft.5|) can be easily performed by using the expression 
for G^ + (q,t — t') given in ( |2.11| ), and by changing the time integration variable to t — t' — 
t" . However, if when computing the above dissipation terms we use the free propagator 
expressions given in Q2.12] ), we would find that they both vanish. We would also obtain the 
wrong results for the temperature dependent vertex and mass corrections coming from the 
second and fourth graphs in ( |2.13| ), as can be explicitly checked. The results would be quite 
different if, instead of free propagators, we use dressed propagators. In fact, consistency 
demands that we correct the propagators for the scalar field by the respective self-energy 
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contributions at least up to two-loops and 0(A 2 ), since this is the order at which we are 
evaluating corrections in perturbation theory. This situation is analogous to the current 
resummation techniques being employed in the improved versions of the electroweak effective 
potential, as can be seen in Ref. 0. We thus write the dressed propagator as 

1 1 



q 2 — m 2 + ie q 2 — m 2 — S(g) + it 
where S(g) is the self-energy contribution, 



(3.6) 



£(?) = O + O + Q + 0(X 3 ) . (3.7) 
In Appendix A we show that the physical propagator Gt + (q,t — t') is then changed to 



e -r{q)\t-t\ 

G^ + (q, t-t')~ [(1 + 2n) cos [u\t - t'\] - i sin [u\t - t'\] 

+ 2f3T(q)n(l + n) sin [uo\t - t'\] + O ( ^ 



(3.8) 



where T(q) is the particle decay width |f26 



. . ImSfo) 

rw = " w (3 ' 9) 

and in ( |3.8|) we used the approximation f3T 1 (see appendix A), which is consistent with 
slow relaxation time-scales. In ( |3.8| ), u = oo(q) and n(u) are now given in terms of the finite 
temperature effetive mass m-r 



r» mT 2 T 2 A 2 T 3 X 2 T 



2 ^ 2 'm 2 



+ ReE(mx) *T m > + A±j - £±- + ^ U ^ j + . . . , (3,0) 

where we have written explicitly the main thermal contributions from each of the terms in 
( p.7[) . The second and third terms in the rhs of ( 3 . 1 0|) are easily obtained. The last term in 



the rhs of ( |3.10| ), associated with the "setting sun" diagram, is explicitly evaluated in p9 

Using the dressed propagator ( p.8|) in the expression for the dissipation term ( |3.4|) and 
performing the integration in t', we obtain, to order A 2 , 
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X \ c (x) [ d 3 x' f dt'<pl(x,t')Im\G+ + ^ ~ 

J J — oo L 



x,x' 



A 2 \ • \ n f d 3 q n(l + n 

—ip c (x,t)ip c (x,t) p 



A 2 , /• d 3 g 1 



(27r)3a; 2 (9)r(g) 
1 + 2n 



+ /?n(l + n) 



+ O ( A 2 — ) . (3.11) 



u> 



The first term in the rhs of ( |3 . 1 1| ) gives the dissipation term, ?7i<£> 2 c , with dissipation 
coefficient r]i given by 

The second term in the rhs in fl3.11|) , clearly gives just the one- loop finite temperature vertex 
correction. In order to obtain (|3.11|) we have performed an expansion to first order in powers 
of T/u, consistent with slowly varying modes. Also, since F ~ 0(X 2 ), we have omitted the 
(9(A 4 ) contributions. The expression for the dissipation coefficient can be further simplified 
if we consider the high temperature limit T ^> mj. As shown in Refs. [[11],^J the high 
temperature limit of T(q) is 

T T> ^ T A2T ' . (3.13) 



Using (|3.13|) in ( 3.12j) , we obtain for rji, in the high temperature limit, 



T»m T 96 / T \ . 
Vi ^ ^ln — , 3.14 
7rT \m T J 

which shows that the dissipation coefficient associated with the multiplicative noise field £i 
is, in this limit, only weakly (logarithmically) dependent on the coupling constant A. 

We can proceed in an analogous way and evaluate Eq. (|3.5|) in order to obtain the 
expression for the dissipation coefficient associated with the second fluctuation (noise) field 



£2, from the first term in the rhs of ( |3.5| ). From the second term we can obtain the two- 
loop mass correction coming from the fourth graph in ( |2.13| ). Substituting Eq. ( |3.8|) for 
G^ + (qj,t—t') in ( |3.5| ) and performing the integration in t', it is possible to show (see appendix 
B) that the dissipation coefficient associated with £ 2 is at least of order X 2 T(qj) ~ 0(X 4 ). 
Therefore, in a weakly interacting model, the dominant contribution to dissipation in the 
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equation of motion for ip c comes from the dissipation term associated with the multiplicative 
noise field, £1. 



B. The effective Langevin-like equation of motion 



Hence, up to two loops and 0(A 2 ), at zero external momentum and within the adiabatic 
approximation we obtain, from (|3.2| ), the following equation of motion for ip c : 



□ + m 2 T (p c (x, t) + ^-ipl(x, t) + 771 t)cp c (x, t) = (p c (x, t) 



(3.15) 



where rji is given by (|3.14|) , uit and At are the renormalized finite temperature mass and 
coupling constant, respectively, obtained from the renormalized effective action, Eq. ( |2.27|) . 
The renormalization of S e e can be defined by the usual introduction of counterterms in the 
initial Lagrangian, Eq. (|2TT|), by writing C — > C+5C, where 5C = \Z{d iii $) 2 — \8m 2 (j) 2 — ^r0 4 , 
with 2, Sm 2 and SX being the wave-function, mass and vertex renormalization counterterms, 
respectively. 5X cancels the logarithmic divergence of the one-loop vertex correction, while 
Z and Sm 2 renormalize the self-energy contribution, Eq. ( |3.7| ). In the high temperature 
limit, in? is given by Eq. ( |3.10|) and is given by 



Ar — A 



3A 



2 r 



T 



+ 



87rm T 8ti 



In + 7' 



+ *(=F) 



(3.16) 



Eq. (|3.15|) can also be written in terms of a finite temperature effective potential Kfr(v 7 c, T), 



96 / T 
□^ C + V^ c ,r) + T^=ln 



(3.17) 



X 2 ttT XTTItJ 

where = d ^\ Vc . 

Note that this equation, apart from the important multiplicative noise source on the rhs, 
is analogous to the one obtained by Hosoya and Sakagami, using quite different methods, 



for the evolution of the thermal average of the scalar field (p c [ITT 



From the equation for the probability distribution for the fluctuation field £1, Eq. 
(2122), we have that the two-point correlation function for £i(x) is given by 
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x 1 



Re 



GV 



(3.18) 



Using ( |2.24|) and ( |3.8|) , we obtain for the two-point correlation function ( |3.18|) the expression 
(at zero external momentum) 

(&(z)£i(z')> = y<J 3 (x-a?) / ( ^)3 4 J (g1 iM") [l + n(u)] + 



+ [1 + 2n(u) + 2n 2 (uj)] cos [2w\t - t'\] + 

+ 2(3T(q)n(uj)[l + n(u)][l + 2n(u)] sin[2^|t - t'\}} e ~ + O (\ 2 ^ 



(3.19) 



which shows that the noise is colored (time dependent), although it is Gaussian distributed. 
Up to order A 2 and for V/u <C 1 , T/T -C 1, we obtain the fluctuation-dissipation relation 

1 



T 



dV<6(*)6(z')>0(*-*') 



(3.20) 



We can also obtain the Markovian limit of fl3.19| ), that is, the limit in which the noise is 
uncorrelated (white). Note that as T — ► oo, T — ► oo, and thus the integrand becomes 
sharply peaked at \t — 1'\ ~ 0. In this limit, we can approximate (|3.19|) by 

d 3 q n(u)[l + n(uj)} 



(2vr)3 2uj\q)Y{q) 
2Tt 1i 5(x - x')5(t - t') , 



(3.21) 



where rji is given by ( |3.12| ). Eq. ( |3.21| ) is the standard expression of the fluctuation- 
dissipation theorem for a Gaussian white noise. 



IV. COUPLING THE SCALAR FIELD TO OTHER FIELDS 

The previous computation of the effective equation of motion for the field configuration 
ifi c can be generalized to include the effects of interactions with other fields. As an example, 
consider the Lagrangian density for the scalar field <fi interacting quadratically with another 
scalar field x, 

m x] = \{d^f + \{d, X f - U[0, X ] , (4.1) 
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with potential V[0, x] given by 

V[*, x] = + ^ + ^X 2 + y + , (4.2) 

where m 2 and fi 2 are positive. This model is a good toy model for several physical cases 
of interest. For example, for some relations among the values of the coupling constants A, 
/ and g 2 (e.g. A ~ 0{g 4 ), f ~ 0{g 2 ) 0), Eq. Q) exhibits the properties of Coleman- 
Weinberg models, for which the quantum corrections coming from integrating out the x field 
break the symmetry in the potential for the scalar field (corrected by the x-loop quantum 
corrections), modifying the original vacuum structure of the model. Also, as pointed out by 
Hu, Paz and Zhang fll3| , (|4~1 ) can mimic, at lowest order (one-loop) in the x-loop quantum 



corrections, a coarse-grained effective model for the scalar field 0, after integrating out the 
X field. In this case, the field would represent the field with components containing the 
long wavelength modes, while x would contain the short wavelength modes, with a cutoff 
determined by some scale A. In inflationary models, would behave as a classical field, 



while x would represent the sub-horizon high frequency modes [[H|,|3^,|nj . The authors in 
|13| thus consider the field x as the quantum bath (at T = 0), allowing them to obtain an 
effective action for the scalar field (the classical action corrected by the x field one-loop 
quantum corrections), where the scalar field is coupled to a noise field, very much like the 
multiplicative noise field £i in Eq. (|2.27|) . Following the results of the last section, the 
generalization of their results to T ^ is relatively simple. Up to one-loop in the x field, 
the effective action T[<p] in ( |2.13|) (also called the influence functional by some authors), will 
be given by 

1 



r[^r^] + -iTr c ln 



D + ^ + gY ■ (4.3) 



Expanding the logarithm in ( |4.3j ) as in (|2.8), up to order g A , we will get expressions analogous 



to the ones in ( |2.15| ) and ( |2.16| ), with (p + , (p_ in the external legs and internal propagators 



for the x field, G™' (x — x') (n, I = +, — ), with expressions just as in (|2.10| ) and ( |2.11| ). By 



changing the field variables <p + , to </?a, V^c as before, the contribution from the x field to 
the effective action for 0, r[(pA,<p c ], Eq. (|2.20|) , will be 
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d xd x' 



d 3 k 1 + 2n A 
(2vr) 3 u; x 



^A(^)<^ c (^)^i(^') + 4<p A (x)<p c (x)<pl(x" 



Im 



-fig d xd x' ip A {x)ip c {x)ip A {x')ip c {x')He 



2 

x,x' 



(4.4) 



where 



2 

x,x' 



is given by an expression analogous to ( 2.24 ). 
The imaginary term in ( |4.4j ), coming from integrating out the \ field (at one-loop order) 
can be rewritten by redefining the fluctuation field £i in ( 2.27Q , such that its probability 
distribution in ( |2.22| ) is changed to 



Pfa] = jVf 1 exp | -~ / d 4 xd 4 x'£ 1 (x) 



—Re 

2 



, + 2g 4 Re 



Ux') ■ (4.5) 



and the two-point correlation function for £i is now given by 



A 2 



-Re 



g: 



+ 2g 4 Re 



G" 



(4.6) 



In the equation of motion for (p c , from ( |4.4p , we will have an additional contribution to the 
dissipation coefficient rji, obtained from a term analogous to (|3.4| ): 

2 



2g 4 / d 4 x'<p 2 c (x')Im 



G: 



, o(t - 1') ~ 

4g A <fi(x,t)(p c (x,t) I dt'(t'-t) 



d 3 q 
(2tt) s 



Im 



G+ + (g,t-r 



G+ + (q,t-t') 



(4.7) 



The first term in the rhs in (|4.7|) gives the contribution to the dissipation coefficient T]i, due 
to the interaction of the scalar field <fi with the x field. The second term in (|4.7|) , together 
with the second term in the rhs in ( f4.4|) , give the corrections of order g 4 and g 2 to the scalar 
field (f) vertex and mass, respectively, due to the x-loop quantum corrections. 

As in ( |3.4| ), in order to obtain a nonvanishing contribution to the dissipation coefficient 
coming from ( |4.7| ), we must consider the dressed propagator G^ + for the \ field, instead of 
the free propagator. G x + (q, t — t') has an expression similar to the one given for the scalar 
field 0, Eq. ([£§) (or Q, for (3T X < 1), where, at 0-{g 4 , f 2 ), the decay width T x can be 
written as, 
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Im£ x (g> x ) 



2w x (g) 

with the imaginary part of the x-held self-energy, from 
of the two two-loop contributions below, 



(4.8) 



given by the imaginary part 



ImSx = Im [ _ _ /_ _\ _ _ ] + Im [ 




X 
/ ~ \ 

4 + 

\ _ / 

X 



(1 



-0qo 



3 r rl 4 k- 



1 + n(kl 



x [g P<p(h)p x ( k 2)p<p(h) + Y^Px( k i)Px( k 2)Px( k z) 



27r) 4 5 4 (g -h-k 2 - k 3 ) , (4.9) 



where p<j,(k) is given by (|A2|) , with m 2 and A corrected by the xToop (T 0) quantum 
corrections. p x (k) is the spectral function for the scalar field Xi with expression analogous 



to the one for the scalar field 0, given by Eq. (|A2|) , but now with T x given by 
given by the solution of u 2 (q) = q 2 + fi 2 + Re£ x (<f, uj x ). 



and u\ 



The high temperature limit of 
field 0, Eq. ( TO ), 



( ->. T»MT 



is analogous to the one for the case with one scalar 



T 2 ( 4 / 2N 
1287rcu x (g) X 12. 



(4.10) 



with = /i 2 + Re£ x (//y). Using these in ( |4.7| ), we obtain an equation of motion for <£> c still 
written as in (|3.15 ), up to two loops and order A 2 in the scalar field <fi and up to one-loopQ 
and order g 4 in the scalar field x- The dissipation coefficient r)\ is given by 



PI 



d 3 q n^l + n^) 
{2'KfumYM 



d 3 q n x (l + n x ) + Q L^_± \ + Q 
(2vr) 3 u\{q)Y x {q) \ u^J \ u Xy 



96 



In 



+ 



\m T J • 12g 4 + f 2ln [fi T/ 



(4.11) 



Up to 2 loops in x the situation would be identical as discussed for the scalar field 0, with results 
similar to the last section and that of appendix B. 
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with the second correction for r]i coming from the interaction in ( |4.1|) . Associated with 
this modified dissipation term there is a modified multiplicative fluctuation (noise) field £1, 
with probability distribution given by ( |4.5| ). For a Coleman- Weinberg potential, we have 
that A ~ Oig 4 ) and / ~ 0(g 2 ), so that the dissipation coefficient is, as in ( |3.14j ), weakly 
dependent on the coupling constants within our approximations. Using the expression for 
the two-point correlation function for £1, Eq. Q4.6Q , both the fluctuation-dissipation relation, 
Eq. (|3.20|) and the Markovian limit expression, Eq. (|3.21|) , still hold. 



V. CONCLUSIONS 

In this work we have studied the nonequilibrium dynamics of a self-coupled scalar field. 
Even though our formalism is in principle applicable in situations far from equilibrium, 
the effective Langevin-like equation we obtained is only adequate to study the approach to 
equilibrium if the initial conditions are not too far from equilibrium. This limitation is essen- 
tially due to the use of perturbation theory and should come as no surprise. However, this 
approach clarifies many important issues concerning nonequilibrium fields and the nature 
of the system-bath coupling. By integrating over fluctuations in order to obtain the effec- 
tive action, it becomes clear in what sense the short wavelength modes can function as the 
thermal bath that drives the longer wavelength modes into equilibrium. In this sense, the 
approximations employed in order to obtain a Langevin-like equation are consistent with this 
system-bath separation; longer wavelength modes have slower dynamics and are responsible 
for the large-distance coherent behavior observed during the approach to equilibrium both 
in the laboratory and in numerical simulations. By going to higher order in perturbation 
theory we were able to obtain the contributions to the noise and dissipation terms coming 
from different diagrams and their relevance to the nonequilibrium dynamics. 

We found that the Langevin-like equation describing the approach to equilibrium both 
for a self-coupled scalar field and for quadratic coupling with other fields is quite different 
from the usual phenomenological form with Gaussian white noise used so far in numerical 
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simulations of the approach to equilibrium in field theory. There are basically three differ- 
ences. The first is that the dominant contribution to the noise is multiplicative; it couples 
quadratically to the field, acting as a "noisy" source to the mass term in the equation of 
motion. The second difference is that even though this noise is still Gaussian distributed, it 
is now non-Markovian; the correlation times depend on the decay width of the fluctuations 
generating the noise. As we show in the text, only in the limit of very high temperatures 
the noise becomes white, as one would naively expect. The final difference has to do with 
the way the dissipation term appears in the equation of motion. Instead of the simple r](f) 
term, we find instead the dissipation "coefficient" depends quadratically on the amplitude 
of the field, t](T)(f) 2 <p. In the high temperature limit for a single scalar field we obtained that 
T](T) ~ (l/T)ln(T/mT), being thus weakly dependent on the coupling constant. This result 
is in agreement with the work of Ref . [|ll[| , which assumed a small departure from equilibrium 



within a kinetic approach. Both results are consistent with linear response theory commonly 
used to obtain transport coefficients in field theories |22| . 



By studying the effects of another scalar field quadratically coupled to <p we were able 
to obtain their different contributions to the noise and dissipation terms in the effective 
equation of motion. Now, the coefficient of the dissipation term depends on ratios of cou- 
plings, as one would expect in more realistic situations, while the noise is still Gaussian 
and multiplicative. In both cases we showed that one can recover a fluctuation-dissipation 
relation. It will be interesting to investigate the implications of this Langevin-like equa- 
tion to the equilibration time-scales during phase transitions, by employing it in numerical 
simulations. Apart from studying the approach to equilibrium from near-equilibirum initial 
conditions, it is possible to use this equation in the study of finite temperature symmetry 
restoration, if one takes into account the effects of expanding about the broken-symmetric 
vacuum. In this connection, it is interesting to note that the coefficient of the dissipation 
term, in the high-temperature limit, displays the typical critical slowing down (poor infra- 
red behavior) observed in many second-order phase transitions; since rj ~ ln(T/m^), as the 
critical temperature is approached from below the temperature corrected mass vanishes and 
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the viscosity diverges logarithmically. We leave as an open question the potential impact 
that a better understanding of nonequilibrium dynamics of field theories will have on our 
current modelling of primordial phase transitions and their possible observational conse- 
quences. However, we believe that interesting physics is lurking behind our present level of 
understanding of nonequilibrium physical processes that took place in the early Universe. 
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APPENDIX A: 

In this Appendix we obtain the expression for the dressed scalar field propagator, Eq. 
The finite temperature, real-time propagator G~f + (q,t — £'), can be written in terms 



(Al) 



of the spectral function p(q,qo) |[26|j22|] , 

G+ + (q,t-t')= f°° ^lp(q,q ){[l + n(q )]9(t-t') + n(q )9(t' - t)} , 
where n(q ) = -p^—; an d the spectral function, for the dressed propagator fl3.6|) , is 

P(g ' go) = ' (q + IT)* - ^{q) ~ (g - if? - uj\q) 

where u(q) is the solution of u: 2 (q) = (f + m 2 + ReS(g, uj) and E(g) is the scalar field self- 
energy, given by (|3.7|) , up to two loops. The spectral function flA2] ) has a peak at go = u(q) 
with a width given by T = T(q), 

ImE(g, a;) 



(A2) 



r(ff) 



Mo) 



(A3) 



For the free propagator, 
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p(q, g ) = i 



(A4) 



(go + ^e) 2 — q 2 — m? (g — ie) 2 — g 2 — m? 

and p(g, go) ^> 27T£:(go)<5(g 2 — m 2 ), where e(g ) = 0(q ) — 8(—q ). Substituting ( [A4] ) in flM|), 
we obtain the free propagator expressions in ( |2.11|) - (|2~T2| ). 

Eq. ( |A2| ) has four poles in the complex go plane: uo ± iT and — c<j ± iT. Using ( |A^ ) in 
(|7n|) and performing the go integration, we obtain 



Gl + (q,t- t') = G>(q,t- t')9{t - H) + G<(q,t- t')6(t' - t) 



(A5) 



where 



G>(q,t- t')= J- {[1 + n(w - ir)] e -*(«-* , X*-*') + n (cu + ir)e^ w+ir ^*-*')} , 



G<{q,t-t')=G>Jq,t'-t) 



(A6) 



The expressions for G, , G^ and G^ + are the same as in ( |2.11| ), but with G^' K given now 
by (|A§). 

T(g) is given in terms of the imaginary part of the self-energy (coming from the third 
graph in Q) by QA3|), where ImS(g) is |2gj22] 



ImS(g) = Im [ 
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-Pqo 



d A ki 



(27r) 4 5 4 (g-A; 1 -A; 2 -A;3) 



(A7) 



The high temperature limit of (|A3|) is given in Refs. |TT| and |29| and we have just quoted 
the final result in the text. 

The expression for G~f + (q,t — t') in ( |A5| ) can also be explicitly written as 

G ™ 1 - ( '> = 2.[cosh^r''cos(/3r)l {smfi ^ C ° sMt " 

+ sin(/?r) sin(w|t - t'|) + z [cos(/?r) - cosh(/?u;)] sin(u;|* - t'|)} . (A8) 



Expanding (|A|) for f3T < 1, we obtain Eq. (gg ). 
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APPENDIX B: 



We estimate here the dissipation coefficient rj 2 associated with the fluctuation field £2, 
obtained from Eq. ( |3.5|) , with G^ + (q, t — t') given by (|X8|) and show that it is subdominant. 
From the first term in the rhs in ( j3.5|) , we get that r] 2 is given by 



02 



A 2 [* ,, f , A f d 3 qi d 3 q 2 



G+ + (qx, t - t')G+ + (q 2 , t - t')G+ + (-qx - q 2 , t - H) 



(Bl) 



From (|A8|) , we can write G~f + (cfj,t — t') as 



Gl + (q 3 ,t-f) = dj + ib 



(B2) 



where aj = a(g},t — t') and bj = b(qj,t — t') are given by the real and imaginary terms of 
Eq. ( |A8| ) , respectively (r^ = T(qj) and ujj = u;(g})) : 
e -v 3 \t-t'\ 



2uij [cosh.{(3ujj) — cos([3Tj)] 



[smh(/3ujj) cos(u;j|t — t'\) + sin(j3Tj) sin(u;,-|t — t'\ 



-r 3 \t-t'\ 



sm(u)j\t — t' |) 



2iUn 



(B3) 



Using (|B2|) in (|Bl|), we get (with q 3 = — qi — q 2 ) 

A 2 /•* ..„., r d 3 gi d\ 2 



A f l ., . r crqi a J q 2 . , , , , , , 

V2 = y J dt(t -t) J 7^317^31 + aia 3 & 2 + a 2 a 3 bx - &1&2O3 



(B4) 



(2tt) 3 (2tt) 3 

Using Eq. ( [B3|) for aj and bj, we can perform the time integration in (|B4|) , by changing the 
time integration variable t' to t — t' = t", and obtain for ([B4]) the expression 

2(1 + 2m)(n 2 - n 3 ) + (1 + 2n 2 )(l + 2n 3 ) - 1 



= a^ r d 3 qi d 3 g 2 1 f £1 + r 2 + r 3 

m 3 J (2tt) 3 (2tt) 3 1 2 

2(1 + 2n 2 )(n 3 - m) + (1 + 2nQ(l + 2n 3 ) - 1 



- u 2 + u 3 y 



+ 



(Ul +LU 2 - LU 3 ) 3 

2(1 + 2n 3 )(n 1 - ra 2 ) + (1 + 2rii)(l + 2n 2 ) - 1 

(— Ux +u 2 + u 3 ) 3 
(1 + 2m)(l + 2n 2 ) + (1 + 2m)(l + 2n 3 ) + (1 + 2w 2 )(l + 2n 3 ) - 1 

(u;i + oo 2 + u 3 ) 3 
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+ /3rim(i + ni) 

+ (n 2 + n 3 ) 



[n 2 - n 3 ) 
1 



(iOi - L0 2 + ^ 3 ) 2 (Ui + L0 2 - 0J 3 ) 2 / 
1 



+ 



+ (3T 2 n 2 (l +n 2 ) 
+ (ni + n 3 



(0Jl + UJ 2 + CJ 3 ) 2 (-CJi + L0 2 + CJ 3 ) 2 / 

1 



+ 



(ni - n 3 ) 
1 



1 



-Ui + CJ 2 + L0 3 ) 2 (Ui + UJ 2 - UJ 3 ) 2 / 

1 



+ 



+ (3T 3 n 3 (l + n 3 ) 
+ (ni + n 2 ) 



(ui +uj 2 + u 3 ) 2 {ui - uj 2 + cj 3 ) 2 / 

1 



+ 



(ni - n 2 ) 
1 



1 



(— + CJ 2 + C^ 3 ) 2 (^1 - UJ 2 + 



y (Ui +UJ 2 + UJ 3 ) 2 {U)\ - Ld 2 + CJ 3 ) 2 y 

The above expression is at least of order A 2 Tj and, since Tj ~ C(A 2 ), we have that r\ 2 ~ 
C(A 4 ). We are thus justified in neglecting its contribution to the effective equation of motion 
to £>(A 2 ). 



(B5) 



27 



REFERENCES 



[1] E. W. Kolb and M. S. Turner, The Early Universe, (Addison- Wesley 1990); 

A. Linde, Particle Physics and Inflationary Cosmology, (Harwood Academic, 1990). 

[2] A. Vilenkin, Phys. Rep. 121, 263 (1985); 

R. H. Brandenberger, Topological Defects and Structure Formation, Lectures given at 
the EPFL, Lausanne, Switzerland, April 1993, Brown report No. Brown-HET-906. 

[3] K. Olive, Phys. Rep. 190, 307 (1990); for natural inflation see, F.C. Adams, J.R. Bond, 
K. Freese, J.A. Frieman, and A. Olinto, Phys. Rev. D47, 426 (1993). 

[4] For a review see, A. G. Cohen, D. B. Kaplan and A. E. Nelson, UC San Diego report 
No. UCSD-PTH-93-02, to appear in Ann. Rev. Nucl. Part. Set., 43. 

[5] G. Boyd, D. Brahm and S.D.H. Hsu, Harvard preprint No. HUTP-92-A027 (1992), in 
press Phys. Rev. D; M. E. Carrington, Phys. Rev. D 45, 2933 (1992); M. Dine, R. 
Leigh, P. Huet, A. Linde, and D. Linde, Phys. Rev. D 46, 550 (1992). J. R. Espinosa, 
M. Quiros, and F. Zwirner, Phys. Lett. B291, 115 (1992); W. Buchmiiller, T. Helbig, 
and D. Walliser, Nucl. Phys. B487, 387 (1993); P. Arnold and O. Espinosa, Phys. Rev. 
D47, 3546 (1993); M. Gleiser and E. W. Kolb, Phys. Rev. D 48, 1560 (1993). 

[6] M. Gleiser and E. W. Kolb, Phys. Rev. Lett. 69, 1304 (1992); G. Gelmini and M. 
Gleiser, Dartmouth report No. DART-HEP-92-06, revised, October 1993; M. Gleiser, 
E. W. Kolb, and R. Watkins, Nucl. Phys. B364, 411 (1991). 

[7] D. Boyanovsky, D. Lee and A. Singh, Phys. Rev. D48, 800 (1993). 

[8] G. F. Mazenko, W. G. Unruh and R. M. Wald, Phys. Rev. D31, 273 (1985). 

[9] A. Albrecht and R. H. Brandenberger, Phys. Rev. D31, 1225 (1985). 

[10] R. Brandenberger, H. Feldman, and J. MacGibbon, Phys. Rev. D 37, 2071 (1988); H. 
Feldman, Phys. Rev. D38, 459 (1988). 

28 



[11] A. Hosoya and M. Sakagami, Phys. Rev. D29, 2228 (1984). 
[12] M. Morikawa, Phys. Rev. D33, 3607 (1986). 

[13] B. L. Hu, J. P. Paz and Y. Zhang, in The Origin of Structure in the Universe, Ed. E. 
Gunzig and P. Nardone (Kluwer Acad. Publ. 1993). 

[14] D. Lee and D. Boyanovsky, Nucl. Phys. B406, 631 (1993). 

[15] D. Boyanosvsky, H.J. de Vega, and R. Holman, University of Pittsburgh report No. 
PITT-93-6, October 1993; J.M. Cornwall and R. Bruinsma, Phys. Rev. D38, 3146 



(1998); B. L. Hu, J. P. Paz and Y. Zhang, in Ref. |13[], and references therein. 



[16] S. Habib, in Stochastic Processes in Astrophysics, Proc. Eighth Annual Workshop in 
Nonlinear Astronomy (1993), Los Alamos report No. LA-UR-93-2105. 

[17] M. Alford, H. Feldman, and M. Gleiser, Phys. Rev. Lett. 68, 1645 (1992); F. J. Alexander 
and S. Habib, Phys. Rev. Lett. 71, 955 (1993); A.I. Bochkarev and Ph. de Forcrand, 
Phys. Rev. Lett. 63, 2337 (1989). 

[18] M. Alford, H. Feldman, and M. Gleiser, Phys. Rev. D47, R2168 (1993). 

[19] M. Alford and M. Gleiser, Phys. Rev. D48, 2838 (1993). 

[20] A. Hosoya, M. Sakagami and M. Takao, Ann. Phys. 154, 229 (1984). 

[21] K. Chou, Z. Su, B. Hao and L. Yu, Phys. Rep. 118, 1 (1985). 

[22] S. Jeon, Phys. Rev. D47, 4586 (1993). 

[23] F.F. Abraham, M. Alford, and M. Gleiser, in progress 

[24] J. Schwinger, Jour. Math. Phys. 2, 407 (1961); L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 
1515 (1964); A. Niemi and G. Semenoff, Ann. Phys. (NY) 152, 105 (1984); Nucl. Phys. 
B230, 181 (1984). 

[25] An accessible introduction to the real-time formalism can be found in R. Rivers, Path 



29 



Integral Methods in Quantum Field Theory, Cambridge University Press, (Cambridge 
1987). 

[26] N. P. Landsman and Ch. G. van Weert, Phys. Rep. 145, 141 (1987). 
[27] P. C. Martin, E. D. Siggia and H. A. Rose, Phys. Rev. A8, 423 (1973). 
[28] M. Morikawa and M. Sasaki, Phys. Lett. 165B, 59 (1985). 
[29] R. R. Parwani, Phys. Rev. D45, 4695 (1992). 
[30] E. J. Weinberg, Phys. Rev. D47, 4614 (1993). 

[31] A. A. Starobinsky, Phys. Lett. 117B, 175 (1982); in Current topics in field theory, 
quantum gravity, and strings, Lecture Notes in Physics, V. 246, Eds. H. J. de Vega and 
N. Sanchez (Springer, Heidelberg 1986). 

[32] S. Rey, Nucl. Phys. B289, 706 (1987). 



30 



